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Abstract 

Frost  heave  modeling  is  presented  os  a  problem  of  small-scale  experimental 
modeling.  Scale  factors  applicable  to  frost  heave  model  testing  are  reviewed, 
and  initial  data  from  frost  heave  experiments  conducted  as  centrifuge  model 
tests  are  presented.  Ongoing  improvements,  modifications  and  future  model 
tests  are  discussed. 
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Initial  Results  From  Small-Scale 
Frost  Heave  Experiments  in  a  Centrifuge 

STEPHEN  A.  KETCHAM  AND  PATRICK  B.  BLACK 


INTRODUCTION 

Frost  heave  is  the  displacement  of  soil 
caused  by  the  formation  of  segregated 
ice  lenses  in  a  freezing  soil  body.  Al¬ 
though  the  process  is  complicated,  a 
number  of  researchers  have  developed 
models  to  describe  and  analyze  the  me¬ 
chanics  of  frost  heave  [see  Kay  and  Per¬ 
fect  (1988)  and  Nakano  (1990)  for  recent 
reviews].  One  of  these  is  Miller's  (1978, 

1980a)  "Rigidice  model." 

In  the  ^gidice  model  the  soil  body 
is  divided  into  three  regions.  The  first  region  com¬ 
prises  already-frozen  soil  that  might  contain  relic 
ice  lenses  but  definitely  contains  a  growing  ice 
lens.  The  second  region,  which  contains  a  grow¬ 
ing  ice  lens  at  one  end  and  unfrozen  soil  at  the 
other,  is  called  the  frozen  fringe  (Miller  1978, 
1980b).  The  third  comprises  unfrozen  soil  and  is 
the  source  of  mass  and  energy  for  the  heaving 
process.  Miller  (1978,  1980b)  calls  this  three-re¬ 
gion  frost  heave  scenario  "secondary  heave."  In 
a  variation  of  secondary  heave,  which  Miller  des¬ 
ignates  "primary  heave,"  there  is  no  frozen  fringe 
but  rather  an  abrupt  boimdary  between  grow¬ 
ing  ice  lens  and  the  unfrozen  soil.  Figure  1  illus¬ 
trates  primary  and  secondary  heaving  conditions. 

Miller  presented  a  mathematical  description 
of  the  Rigidice  model  in  1978.  He  followed  with 
finite-element  (O'Neill  and  Miller  1985)  and  fi¬ 
nite-difference  implementations  (Black  and  Miller 
1985),  comparisons  to  experimental  data  (Black 
and  Miller  1985)  and  a  scaling  analysis  of  the  gov¬ 
erning  equations  of  the  model  (Black  1985,  Miller 
1990).  Of  particular  interest  for  this  paper  is 
Miller's  scaling  analysis.  This  analysis  and  the 


a.  Primary  heave.  b.  Secondary  heave. 

Figure  1.  Frost  heave  scenarios. 

eling  of  frost  heave  would  best  be  conducted  as 
a  body  force  analog  and  that  a  small-scale  frost 
heave  test  conducted  in  a  centrifuge  would  be 
appropriate  for  this  purpose. 

An  objective  of  our  ongoing  work  is  to  experi¬ 
mentally  confirm  Miller's  conclusion  and  the  scale 
factors  inherent  in  his  scaling  analysis.  This  re¬ 
port  presents  a  brief  overview  of  frost  heave  cen¬ 
trifuge  modeling  from  two  perspectives:  conven¬ 
tional  geotechnical  centrifuge  modeling  and  scal¬ 
ing  of  the  differential  equations  of  the  Rigidice 
model.  We  describe  two  initial  model  experiments 
in  which  silt  specimens  were  partially  frozen  in 
8-^  and  32-^  centrifugal  fields.  The  experiments 
were  performed  to  provide  a  basis  for  equipment 
development  and  further  experiments.  Plans  for 
the  equipment  and  future  tests  are  presented  as 
well. 


CENTRIFUGE  MODELING 
BACKGROUND 

Physical  modeling  is  the  testing  of  a  small-scale 


distinct  role  of  gravitational  body  force  in  the  heav-  model  so  that  the  response  of  a  prototype  can  be 

ing  process  led  Miller  to  conclude  that  scale  mod-  inferred  from  the  model  response.  A  model  is 


\ 


Table  1.  Scale  factors  for  geotechnical  centrifuge  mod¬ 
eling.  (After  Scott  and  Morgan  1977,  Croce  et  al.  1985, 
Savidou  1988.) 

Quantity  Prototype  Model 


Linear  dimension  and  displacement  1  1/N 

Area  1  l/N^ 

Volume  1  1/N^ 

Mass  density  1  1 

Mass  1  1/N^ 

Acceleration  1  N 

Stress  1  1 

Force  1  1/N^ 

Strain  1  1 

Energy  density  1  1 

Energy  1  1/N^ 

Temperature  1  1 

Time  (for  viscous  force  similarity)  1  1 

Time  (for  inertial  force  similarity)  1  1/N 

Time  (for  seepage  force  similarity)  1  1/N^ 


said  to  be  similar  to  the  prototype  when  each  sig¬ 
nificant  engineering  variable  of  the  model  is  re¬ 
lated  by  a  proportionality  or  scale  factor  to  the 
corresponding  variable  of  the  prototype.  Scale 
factors  are  governed  by  the  physics  and  model 
laws  of  the  problem.  They  are  used  to  design  the 
model  and  to  interpret  the  measured  model  re¬ 
sponse  as  the  prototype  response. 

Physical  models  of  civil  engineering  structures 
are  often  designed  so  that  the  model  will  be  geo¬ 
metrically  similar  to  the  prototype,  so  that  sig¬ 
nificant  forces  in  the  model  are  proportional  to 
the  forces  in  the  prototype,  and  so  that  the  model's 
response  to  loads  will  be  kinematically  similar 
to  that  of  the  prototype.  Geometric  similarity 
means  that  all  parts  of  the  model  have  the  same 
shapes  as  the  corresponding  parts  of  the  proto¬ 
type,  and  kinematic  similarity  means  that  the 
motions  of  the  model  are  similar  to  the  motions 
of  the  prototype  at  corresponding  times.  When  all 
net  forces  are  proportional,  dynamic  similarity  is 
said  to  exist  (Langhaar  1951).  In  general,  the  me¬ 
chanical  response  of  a  structural  model  to  forces 
is  observed  and  interpreted  as  the  prototype  re¬ 
sponse. 

Centrifuge  modeling  is  a  physical  modeling 
technique  in  which  the  weight  stresses  of  a  struc¬ 
ture  are  simulated  by  the  placement  of  a  small- 
scale  model  in  a  centrifugal  field.  This  technique 
is  of  proven  benefit  for  modeling  soil  structures, 
because  the  form  and  magnitude  of  the  soil  re¬ 
sponse  are  often  greatly  dependent  on  weight¬ 
generated  effective  stresses  and  because  the  load¬ 
ing  can  be  dominated  by  weight  loads.  As  de¬ 
scribed  in  Appendix  A,  when  the  prototype  soil 


is  used  in  the  model,  model  laws  of  soil  mechan¬ 
ics,  together  with  requirements  for  geometric  and 
kinematic  similarity,  lead  to  scale  factors  for  cen¬ 
trifuge  modeling.  Several  are  shown  in  Table  1  as 
ratios  of  the  model  quantity  to  the  prototype  quan¬ 
tity  for  a  model  that  is  1/N  the  size  of  the  proto¬ 
type.  Following  Croce  et  al.  (1985),  a  number  of 
these  are  derived  in  Appendix  A. 

Many  of  the  scale  factors  listed  in  Table  1  are 
basic  structural  modeling  scale  factors  and  can 
be  derived  from  the  mechanical  considerations 
(Langhaar  1951).  The  scale  factor  for  acceleration, 
which  is  N,  can  be  derived  from  the  consider¬ 
ation  that  weight  forces  should  be  scaled  as  other 
forces,  and  it  implies  the  need  for  an  increased 
gravity  field  (Croce  et  al.  1985).  The  centrifuge  is 
used  to  approximate  this  condition  by  subject¬ 
ing  a  model  to  the  constant  angular  velocity  Q. 
according  to 


where  g  ~  acceleration  due  to  Earth's  gravity 

r  =  representative  radius  of  the  model 
from  the  axis  of  rotation 
rQ?-  =  radial  acceleration  of  the  model  at  this 
radius. 

In  this  context,  N  can  be  considered  as  a  nominal 
gravity  level: 


and  the  model  can  be  considered  to  be  subjected 
to  an  inertial  field  equivalent  to  N  gravities,  or  N 
g.  As  described  by  Langhaar  (1951),  if  self- weight 
stresses  have  a  negligible  influence  on  the  proto¬ 
type  response,  or  if  partial  dynamic  similarity  is 
acceptable  to  the  modeler,  the  principles  discussed 
here  for  centrifuge  modeling  and  the  scaling  re¬ 
lations  of  Table  1  can  be  applied  to  models  tested 
at  1  g.  However,  because  self-weight  stresses  typi¬ 
cally  cannot  be  neglected  for  geotechnical  proto¬ 
types,  centrifuge  modeling  has  become  a  conven¬ 
tional  technique. 

As  indicated  in  Table  1,  there  are  different  time 
scales  for  the  forces  of  viscous,  inertial  and  seep¬ 
age  phenomena.  As  a  result,  time  scale  conflicts 
can  occur  for  certain  modeling  problems,  and  dy¬ 
namic  similarity  can  be  impossible  to  achieve. 
Thus  the  experimenter  must  consider  the  limita¬ 
tions  imposed  by  the  model  laws  and  scale  fac¬ 
tors  when  designing  a  model  experiment.  Appen¬ 
dix  A  presents  derivations  of  these  time  scales. 
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The  use  of  prototype  soils  in  centrifuge  mod¬ 
els  is  an  accepted  practice  and  is  an  attempt  to 
achieve  structural  response  similarity  of  the  model 
and  prototype  by  ensuring  that  the  constitutive 
response  of  the  model  will  be  as  close  as  possible 
to  the  prototype  material  response.  Langhaar 
(1951)  and  Schmidt  and  Holsapple  (1980)  have 
suggested  that  when  a  model  is  made  of  the  ma¬ 
terial  of  the  prototype,  when  there  are  no  time 
scale  conflicts,  and  when  the  material  response 
is  independent  of  model  size,  prototypes  with 
general  constitutive  behavior  such  as  nonlinear 
elasticity,  plasticity  and  fracture  can  be  modeled 
by  following  the  scaling  relations  of  Table  1.  The 
use  of  prototype  soils  within  typical  scale  mod¬ 
els  (e.g.,  1/25  to  1/100  scale)  does  not  generally 
result  in  scale  effects  imless  soils  with  particles 
larger  than  sands  are  tested.  It  is  common  in  cen¬ 
trifuge  testing  to  ensure  that  the  parts  of  models 
are  much  larger  than  the  soil  particle  sizes  so  that 
particle  size  does  not  result  in  adverse  scale  ef¬ 
fects.  It  is  also  common  to  test  different  scale 
models  of  the  same  prototype  to  ensure  that  scale 
effects  are  insignificant. 

To  correctly  model  heat  transfer  processes  in 
saturated  ground,  the  model  and  prototype  must 
be  thermally  similar,  in  addition  to  being  geometri¬ 
cally,  kinematically  and  dynamically  similar.  Ther¬ 
mal  similarity  includes  the  similarity  of  internal 
energy  and  energy  density,  as  evidenced  by  similar 
temperature  changes  and  differences  at  homolo¬ 
gous  points.  For  centrifuge  modeling,  Savidou 
(1988)  has  shown  that  the  scale  factor  for  tem¬ 
perature  change  is  1  and  that  the  time  scale  fac¬ 
tor  for  conductive  and  convective  heat  transfer 
through  a  soil  body  is  1  /N^,  which  is  the  same  as 
the  time  scale  factor  for  seepage  processes. 

FROST  HEAVE  MODELING 

When  the  differential  equations  of  a  particu¬ 
lar  mechanical  problem  or  physical  process  are 
known  to  the  modeler,  a  scaling  analysis  can  be 
performed  to  determine  techniques  for  appropri¬ 
ate  small-scale  modeling  of  the  problem.  Scaling 
is  a  method  of  dissolving  passive  physical  pa¬ 
rameters  into  active  variables  in  a  differential 
equation  (Miller  1980a).  Any  solution  to  the  scaled 
equation  is  also  a  solution  to  any  other  distinct 
physical  system  that  can  obtain  the  same  scaled 
variables  and  scaled  boundary  conditions  from 
its  distinct  physical  parameters.  The  earliest  scaling 
analysis  of  soils  containing  water  and  its  vapor 


was  conducted  by  Miller  and  Miller  (1956).  They 
presented  a  set  of  scaled  relationships  describ¬ 
ing  the  isothermal  flow  of  liquids  in  the  unsatur¬ 
ated  soil,  and  they  concluded  from  length,  stress 
and  time  scale  factors  that  a  centrifuge  experi¬ 
ment  would  be  appropriate  for  small-scale  mod¬ 
eling  of  this  flow  (Miller  and  Miller  1955).  Miller's 
1990  scaling  analysis  of  the  Rigidice  equations 
was  an  extension  of  these  earlier  investigations 
of  scaling  laws  for  unsaturated  porous  flow.  An 
expansion  of  this  analysis  is  presented  in  Ap¬ 
pendix  B. 

Like  the  previous  scaling  investigations. 
Miller's  Rigidice  analysis  was  based  on  scaling 
the  effects  of  surface  tension  and  viscous  flow, 
but  it  also  incorporated  heat  transfer,  phase  change 
of  the  pore  water,  and  ice  lens  initiation  and 
growth.  The  analysis  yielded  reduced  equations 
that  are  consistent  with  the  techniques  of  con¬ 
ventional  centrifuge  modeling  and  the  scale  fac¬ 
tors  in  Table  1.  It  provided  a  complete  theoreti¬ 
cal  basis  for  small-scale  frost  heave  modeling  in 
a  geotechnical  centrifuge.  Miller  (1990)  empha¬ 
sized  that  the  ratio  of  model  time  to  prototype 
time  for  the  freezing  process  in  a  centrifuge  model 
would  be  1/N^,  and  he  cited  Pokrovsky  and 
Fyodorov  (1969)  as  the  first  to  suggest  the  use  of 
a  geotechnical  centrifuge  for  modeling  frozen  and 
freezing  ground  structures.  On  the  basis  of  Miller 's 
earlier  unpublished  discussions  on  model  laws 
for  soil  freezing.  Black  (1985)  presented  a  set  of 
scaling  relationships  for  frost  heave  centrifuge 
modeling,  recommended  a  simple  model  study 
of  frost  penetration  and  heave,  and  demonstrated 
that  a  great  time  savings  would  be  obtained  over 
a  full-scale  experiment. 

EXPERIMENTAL  TECHNIQUE 

Our  preliminary  soil  freezing  experiments  were 
conducted  on  an  International  Equipment  Com¬ 
pany  Model  PR-2  centrifuge.  The  device  is  a  small 
refrigerated  centrifuge  manufactured  for  appli¬ 
cations  in  biology,  bacteriology,  physiology  and 
biochemistry.  It  has  a  nominal  operating  radius 
of  0.2  m  for  the  rotor  and  cups  used  in  our  test 
configuration.  During  the  experiments  the  centri¬ 
fuge  was  located  in  a  coldroom  with  an  ambient 
temperature  of  approximately  4.5°C. 

Two  insulated  sample  containers  were  made 
to  allow  simultaneous  testing  of  saturated  soil 
samples.  In  each  container  a  Plexiglas  cylinder 
enclosed  the  soil  sample.  The  cylinder  was  moimt- 
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ed  on  a  porous  plastic  base  to  allow  water  flow  to 
the  sample  from  an  exterior  annular  reservoir.  Ex¬ 
cept  for  the  upper  end  of  the  Plexiglas  cylinder, 
the  cylinder  and  reservoir  were  surrounded  by 
insulation.  The  upper  end  of  the  cylinder  was  left 
open  to  the  air.  The  insulated  sample  containers 
were  made  to  fit  into  cylindrical  cups  that  were 
supplied  with  the  centrifuge.  These  cups  have 
inside  dimensions  of  98  mm  in  diameter  and  120 
mm  in  height. 

The  test  samples  comprised  a  saturated  New 
Hampshire  silt  compacted  in  the  Plexiglas  cylin- 


Figure  2,  Sample  containers  and  sample  preloading  equip¬ 
ment. 


der  of  the  assembled  sample  container.  The  nomi¬ 
nal  sample  dimensions  were  10  mm  in  height  and 
34  mm  in  diameter.  We  prepared  each  sample  at 
room  temperature  by  mixing  12.5  g  of  the  New 
Hampshire  silt  with  water  to  form  a  slurry  with 
an  80%  water  content.  The  slurry  was  poured  into 
the  Plexiglas  cylinder,  preloaded  and  consolidated 
with  a  vertical  stress  of  14.1  kPa.  Free  water  ex¬ 
pelled  at  the  sample  surface  was  removed.  The 
dry  mass  density  of  the  consolidated  samples  was 
approximately  2040  kg/m^.  The  sample  contain¬ 
ers  and  sample  preloading  equipment  are  shown 
in  Figure  2. 

Following  the  sample  preparation,  the  in¬ 
sulated  sample  containers  were  taken  to  the 
coldroom,  placed  in  the  centrifuge  cups,  and 
immediately  tested.  Figure  3  shows  the  sample 
containers  in  the  centrifuge.  The  radius  from 
the  centrifuge  axis  to  the  sample  bottom  was 
0.196  m.  To  achieve  nominal  gravity  levels  of 
8  and  32  g  at  desired  reference  radii  in  the  tests 
(Cooke  1991),  the  centrifuge  was  operated  at 
rotational  velocities  close  to  193  and  387  rpm, 
respectively.  The  effect  of  Earth's  gravity  on  the 
equilibrium  orientation  of  the  centrifuge  cups  was 
considered  in  our  calculation  of  these  velocities. 
The  experiments  were  conducted  with  mini¬ 
mal  control  of  the  temperature  of  the  sample 
and  with  minimal  measurement  of  the  sample 
response.  Our  hope  for  these  preliminary  ex¬ 
periments  was  simply  that  we  could  freeze  the 
samples  from  the  top  down.  We  operated  the 
centrifuge  with  the  centrifuge  chamber  closed 
to  the  coldroom  atmosphere,  and  we  main¬ 
tained  the  chamber  temperature  close  to  our 
desired  temperature  of  “2°C  by  manually  turn¬ 
ing  on  and  off  the  refrigeration  system  of  the 


Figure  3.  Sample  containers  in  the  centrifuge. 
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centrifuge.  Manual  operation  was  necessary  be¬ 
cause  of  the  crude  response  of  the  centrifuge's 
refrigeration  control  system.  The  duration  of  the 
experiments  was  approximately  five  hours.  A  ther¬ 
mocouple  located  at  the  top  of  the  closed  centri¬ 
fuge  chamber  provided  the  only  temperature 
measurement  for  the  tests.  There  was  no  measure¬ 
ment  of  the  temperature  of  the  sample.  In  addi¬ 
tion  there  was  no  measurement  of  the  sample 
heave  as  a  fimction  of  time. 


RESULTS 

The  experimental  data  consist  of  measurements 
of  the  temperature  in  the  centrifuge  chamber,  mea¬ 
surements  of  the  height  of  the  samples  at  the  be¬ 
ginning  and  end  of  the  tests,  and  observations 
and  photographs  of  the  samples  following  the 
tests.  Table  2  presents  the  measurements  for  the 


8-^  and  32-^  tests  along  with  the  corresponding 
1-g  prototype  values  calculated  using  either  the 
Table  1  scale  factors  or  the  scaling  process  de¬ 
scribed  in  Appendix  B.  The  temperature  records 
reveal  that  the  chamber  temperature  varied  by 
±0.25°C  from  the  desired  temperature  of  -2°C. 
Figures  4  and  5  show  the  8-g  and  32-^  samples, 
respectively,  before  and  after  their  extrusion  from 
the  cylinder  and  after  crudely  cutting  them  to 
examine  the  interior  ice  lens  configuration.  The 
air  temperature  surrounding  the  samples  during 
cutting  was  4.5°C.  As  a  result  the  figures  show 
the  samples  in  progressive  stages  of  melting. 

In  the  five-hour  test  period,  only  the  upper  5- 
8  mm  of  the  8-^  and  32-^  samples  froze.  This  was 
a  desired  result,  as  it  confirms  that  the  freezing 
occurred  from  the  top  down.  This  was  important 
since  we  had  no  other  means  of  monitoring  the 
progress  of  the  freezing  front.  That  the  8-^  and 
32-^  samples  would  freeze  to  about  the  same  depth 


b.  After  extrusion  from  the  cylinder. 


c.  After  cutting  to  examine  the  interior  ice  lens 
configuration. 

Figure  4.  Sample  used  in  the  8-g  experiment. 
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Table  2.  Centrifuge  test  conditions  along  with  measured  model  and  calculated 
prototype  results  for  New  Hampshire  silt. 


Quantity 

1/8-scale  experiment 

1/32-scale  experiment 

Model 

Prototype 

Model 

Prototypi 

Acceleration  (g) 

8 

1 

32 

1 

Dry  bulk  density  (kg/m^) 

2040 

2040 

2040 

2040 

Sample  height  (mm) 

10 

80* 

10 

320* 

Sample  diameter  (mm) 

34 

272* 

34 

1088* 

Upper  surface  temperature  (°C) 

-2 

-2+ 

-2 

-2+ 

Duration  of  test  (hr) 

5 

320** 

5 

5120** 

Frost  penetration  (mm) 

5-8 

40-64* 

5-8 

160-256* 

Heave  (mm) 

1-2 

8-16* 

3 

96* 

*Table  1  or  eq  B18 
^Table  1  or  eq  BIO 
**Table  1  or  eq  B24  with  B18 


c.  After  cutting  to  examine  the  interior  ice  lens  con¬ 
figuration. 

Figure  5.  Sample  used  in  32-g  experiment. 
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was  also  an  expected  result,  since  an  inertial  field 
should  not  directly  affect  the  heat  transfer.  From 
the  sample  height  measurements  we  estimated 
the  average  heave  of  each  of  the  two  S-g  samples 
to  be  approximately  1-2  mm  and  the  average 
heave  of  each  of  the  two  32-g  samples  to  be  ap¬ 
proximately  3  mm. 

The  S-g  test  samples  froze  imder  a  stress  field 
that  is  approximately  similar  to  the  stress  field  of 
an  8-cm  column  of  sod  at  1  g,  and  the  32-^  test  samples 
froze  imder  a  stress  field  approximately  similar  to 
a  32-cm  column  of  soil  at  1  g.  The  differences  in  the 
configurations  of  the  interior  lenses  of  the  8-^samples 
and  die  32-g  samples  were  the  most  interesting  re¬ 
sult  of  the  experiments.  In  the  S-g  samples  the  ice 
lenses  did  not  form  in  distinct  horizontal  layers. 
The  ice  clearly  formed  thin  (<  0.5  mm)  lenses  that 
separated  the  soil  particles,  but  the  lenses  appeared 
randomly  oriented.  In  the  32-^  samples,  however, 
a  horizontal  bias  to  the  lens  orientation  was  evi¬ 
dent.  These  observations  are  consistent  with  an 
explanation  that  might  be  used  for  similar  observa¬ 
tions  in  the  upper  4  cm  and  the  upper  16  cm  of  a 
silt  column  in  1-g  freezing  tests:  that  a  randomly 
oriented  lens  configuration  can  develop  in  shal¬ 
low  soils  that  do  not  have  a  significant  vertical 
overburden  stress,  and  a  horizontally  oriented  lens 
pattern  can  develop  in  deeper  soils  that  have  a 
greater  overburden  stress. 


The  results  of  the  experiments  are  crude,  yet 
they  provide  two  pieces  of  information  for  the 
design  of  equipment  for  further  experiments.  First, 
the  indication  that  the  freezing  occurred  from  the 
top  down  in  the  samples  suggests  that  we  can 
proceed  with  a  one-dimensional  cylindrical  sam¬ 
ple  configuration  in  our  small  centrifuge  tests, 
and  second,  the  apparent  consistency  of  the  lens 
configurations  of  the  S-g  and  32-g  samples  with 
corresponding  prototype  behavior  suggests  that 
with  control  of  the  experiments  the  applicability 
of  the  centrifuge  for  frost  heave  modeling  can  be 
investigated. 

CONCLUSIONS  AND 
FUTURE  EXPERIMENTS 

A  design  of  a  simple  frost  heave  cell  for  the 
small  centrifuge  is  illustrated  schematically  in 
Figure  6.  The  cell  will  be  instrumented  with  lower 
and  upper  heaters  and  with  temperature  and  dis¬ 
placement  measurement  sensors.  The  lower  heater 
will  be  a  thin  heating  pad  mounted  onto  the  bot¬ 
tom  of  the  centrifuge  cup,  and  the  upper  heater 
will  be  a  cartridge  heater  mounted  in  a  radial-fin 
heat  sink.  A  fan  will  be  installed  above  the  heater 
to  circulate  air  past  the  heat  sink  fins  to  the  up¬ 
per  sample  surface.  The  upper  heater  will  not  be 


Electric  Fan 


Rotation  Axis 
of  Centrifuge  Cup 


% 


\wm 
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Radial 

Fin 

Heat 
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Disk 


Figure  6.  Design  of  a  simple  frost  heave  cell  for  the  small  centrifuge.  The 
diagram  show  the  cross  section  of  the  inside  of  the  centrifuge  chamber  during 
operation;  the  centrifuge  rotor  arm  is  not  shown. 
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in  contact  with  the  sample  and  thus  will  not  load 
the  sample.  Both  heaters  will  be  controlled  using 
feedback  controllers.  An  electrical  slip  ring  as¬ 
sembly  will  be  installed  on  the  centrifuge  axis  to 
allow  power  and  instrumentation  signals  to  be 
passed  to  and  from  the  cell  while  the  centrifuge 
is  operating.  The  centrifuge  will  be  operated  in  a 
coldroom  at  a  temperature  below  freezing  (e.g., 
-5°C),  and  the  heaters  will  be  controlled  to  main¬ 
tain  upper  and  lower  sample  surfaces  at  desired 
temperatures.  Thin  sectioning  and  image  analyz¬ 
ing  techniques  will  be  used  to  examine  the  lens 
ice  distribution,  orientation  and  thickness. 

The  modeling  at  different  scales  or  accelera¬ 
tion  fields  will  be  used  to  examine  the  validity  of 
our  frost  heave  modeling.  Specifically  a  series  of 
different  scale  models  of  the  same  prototype  will 
be  tested,  and  the  model  displacement  vs.  time 
records  will  be  interpreted  to  see  if  each  of  the 
prototype  heave  vs.  time  predictions  are  essen¬ 
tially  the  same.  If  the  one-dimensional  heave  pro¬ 
cess  can  be  modeled,  then  other  tests  will  be  con¬ 
ducted  regarding  heave  forces  and  the  interaction 
of  heaving  soil  with  foundation  and  other  struc¬ 
tural  members. 

Additional  analysis  will  compare  experimen¬ 
tal  measurements  to  predictions  using  a  recent 
version  of  the  Rigidice  model  (Black  1994).  While 
the  exact  location  and  thickness  of  individual  ice 
lenses  in  the  model  and  prototype  will  most  likely 
not  be  predictable,  the  total  amount  of  heave,  pen¬ 
etration  and  fluxes  should.  This  will  be  the  topic 
of  future  discussions. 
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APPENDIX  A:  SOME  SCALE  FACTORS  FOR 
CONVENTIONAL  CENTRIFUGE  MODELING 


PHYSICAL  MODELING  CONCEPTS 

Following  the  presentation  of  Langhaar  (1951), 
consider  a  physical  variable  of  a  prototype /p(Xp, 
j/p,  Zp,  fp)  and  the  corresponding  model  quantity 
Vm'  ^m'  ^m)'  where  X,  y  and  z  are  coordi¬ 
nate  measures  of  a  point  within  the  structure,  f  is 
time,  and  the  subscripts  p  and  m  refer  to  the  pro¬ 
totype  and  model,  respectively.  The  function  /„ 
is  said  to  be  similar  to  the  function /p  if  the  ratio 
/m//p  is  constant  for  homologous  points  and  ho¬ 
mologous  times.  The  constant  ratio  /  =  /„,/ L  is 
called  the  scale  factor  for  the  function /.  Scale  fac¬ 
tors  can  be  derived  from  model  laws,  i.e.  ratios 
of  physical  laws  for  the  model  and  prototype. 
When  similarity  is  achieved,  a  model's  response 
is  interpreted  using  the  scale  factors  to  infer  the 
response  of  the  corresponding  prototype. 

The  ideal  for  physical  modeling  is  complete 
similarity,  which  includes  geometric,  kinematic, 
thermal  and  dynamic  similarity.  Complete  simi¬ 
larity  is  usually  not  achieved  for  a  model  test. 
Scale  factor  conflicts  and  scale  effects,  i.e.,  dis¬ 
turbing  influences  that  are  associated  with  the 
small  scale  of  the  model,  can  limit  similarity  to  a 
partial  similarity. 

When  a  model  is  made  of  the  material  of  the 
prototype,  and  when  there  are  no  scale  factor  con¬ 
flicts  or  disturbing  scale  effects,  then  the  model 
material  will  have  the  same  constitutive  behav¬ 
ior  as  the  prototype  material,  and  prototypes  with 
general  constitutive  behavior  such  as  nonlinear 
elasticity  and  plasticity  can  be  modeled. 


CROCE'S  DERIVATION  OF 
SCALE  FACTORS  FOR  CONVENTIONAL 
GEOTECHNICAL  MODELING  USING 
SCALAR  MODEL  LAWS 


scalar  model  laws.  Closely  following  the  analy¬ 
sis  of  Croce  et  al.  (1985),  this  approach  is  presented 
here.  Two  assumptions  are  made:  that  soil  can 
be  treated  as  a  continuum  and  that  soil  proper¬ 
ties  are  not  affected  by  a  change  in  acceleration. 

D)mamic  similarity  requires  that  all  forces  (and 
all  kinds  of  forces)  have  the  same  scale  factor.  The 
forces  of  interest  include  the  weight  force  F^,  the 
external  force  Fg,  the  viscous  force  Fy,  the  inertia 
force  Fj  and  the  seepage  force  Fg. 

Scale  factors  for  length  I,  mass  m  and  time  t 
are 
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(Al) 

(A2) 

njp 

Kt--p. 

(A3) 

P 


Thus  the  scale  factors  for  area  A  and  volume  V 
are 

=  (A4) 

Ap 

and 

Kv=Kf=^.  (A5) 

Since  the  same  material  is  used  in  the  model  and 
prototype  at  the  same  mass  density,  the  scale  factor 
for  mass  is  equal  to  the  scale  factor  for  volume: 

^rxv=^\=—- 

Wp 

For  identical  effective  stresses  in  the  model  and 
prototype. 


The  objective  of  conventional  small-scale 
geotechnical  modeling  is  to  achieve  a  model  re¬ 
sponse  that  is  similar  to  the  mechanical  behavior 
of  the  prototype,  and  the  approach  is  to  use  the 
prototype  soil  in  a  geometrically  similar  model. 
Thus  it  is  important  to  ensure  that  the  effective 
stresses  in  the  model  are  the  same  as  in  the  pro¬ 
totype  at  homologous  points  and  times.  The  scale 
factors  for  these  conditions  can  be  derived  using 


u 


^  (A7) 

C^P  CJp  Wp  ' 


where  a',  a  and  u  are,  respectively,  the  effective 
stress,  the  total  stress  and  the  excess  pore  water 
pressure,  including  the  capillary  pressure  caused 
by  surface  tension.  Considering  eq  A4  and  A7  and 
the  definition  of  stress,  it  follows  that  the  scale 
factor  for  all  forces  f  should  be 
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JCp=K2=^.  (A8) 

From  the  definition  of  weight  force,  the  weight 
force  model  law  and  scale  factor  Kv  is 

^Fw=T^  = - ^  =  ^1  (A9) 


Pi  =  ma  (A16) 

where  a  is  the  acceleration.  The  model  law  and 
scale  factor  JCp.  for  the  inertia  force  is 


r'i  nt  Cl 

V  _ 

- 

^ip  “p 


(A17) 


From  the  definition  of  acceleration  it  follows  that 


where  flg  is  the  acceleration  due  to  gravity.  Con¬ 
sidering  eq  A6,  A9  and  the  relation  iCp  =Kj;  =Kf 
from  eq  A8, 


%L  =  1. 
«gn  ^1 


(AlO) 


Equation  AlO  implies  that  geotechnical  model¬ 
ing  requires  a  technique  that  elevates  the  gravi¬ 
tational  acceleration  of  the  model. 

The  scale  factor  for  an  external  force  must 
satisfy  eq  A8,  and  only  technical  constraints  should 
limit  this. 

A  viscous  force  acting  on  a  small  area  A  can  be 
defined  as 


Kf  ap  ■ 

From  eq  A6  and  A18,  eq  A17  becomes 

Because  JCp.  must  equal  iCp  =  ,  then 


(A19) 


(A20) 


That  is,  the  time  scale  factor  is  equal  to  the  length 
scale  factor  if  the  inertia  force  is  to  be  scaled  like 
the  weight  force. 

A  seepage  force  can  be  written  as 


Py  ~  Ps 


F,  =  zWp 


where  [ig  is  the  viscosity  of  the  soil  skeleton  and 
di?/dn  is  the  velocity  gradient.  The  model  law 
and  scale  factor  for  the  viscous  force  Kp  can  then 
be  written 

Pv  M'S  dv  A 


^3  dUp  dn 


(A12) 

m  -^p 


From  the  definition  of  velocity  it  follows  that 

=  (A13) 

From  eq  Al,  A4  and  A13  and  the  fact  that  viscos¬ 
ity  is  independent  of  gravity,  Kp^  can  be  written 

Kv  =^.  (AU) 


Because  Kp  must  equal  Kp  =  Kf,  then 


K^  =  l 


(A15) 


That  is,  the  time  scale  factor  is  1  if  the  viscous 
force  is  to  be  scaled  like  the  weight  force. 

An  inertia  force  can  be  written  as 


where  Wp  is  the  weight  of  the  fluid  phase  and  i  is 
the  hydraulic  gradient  defined  by 

i  =  j  (A22) 

in  which  v  is  the  superficial  velocity  of  the  fluid 
and  k  is  Darcy's  coefficient  of  permeability.  From 
eq  A21  and  A22  the  scale  factor  Xp^  for  the  seep¬ 
age  force  can  be  written  as 

Ps  V  K 

Xp  =  (A23) 

■^s  p  77  Ir  W  '  ' 

^Sp  ’^^Fp 

Equation  A8  implies  that 


-la-  =  Kp  = 


(A24) 


Then,  from  eq  A13  and  A24,  eq  A23  becomes 


(A25) 


A  general  expression  for  the  coefficient  of  per¬ 
meability  is  (Mitchell  1976) 
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where  and  are,  respectively,  the  density 
and  the  viscosity  of  the  fluid  phase,  and  k  is  de¬ 
fined  as  the  absolute  permeability,  which  depends 
only  on  the  geometry  of  the  soil  skeleton.  Because 
Pw,  Mw  ^  depend  on  gravity,  it  follows 

that 


_  8n 


(A27) 


Thus  eq  A25  becomes 


V  —  I 

-T"- 


(A28) 


Because  Kp^  must  equal  Kp  =  K^,  then 
K,=K^. 


That  is,  the  time  scale  factor  is  equal  to  the  length 


scale  factor  squared  if  the  seepage  force  is  to  be 
scaled  like  the  weight  force. 

CENTRIFUGE  TECHNIQUE 

The  centrifuge  is  used  to  approximate  the  con¬ 
dition  of  eq  AlO  by  subjecting  a  model  to  the  con¬ 
stant  angular  velocity  Q  such  that  the  following 
approximate  relation  holds: 

fSnL  =  J_  =  l£^  (A30) 

^gp  "gp 

where  ag^  is  the  acceleration  due  to  Earth's  grav- 
ity,  r  is  a  representative  radius  of  the  model  from 
the  axis  of  rotation  (Cooke  1991),  and  rQ^  is  the 
radial  centripetal  acceleration  of  the  model  at  this 
radius.  In  the  centrifuge  modeling  literature,  N 
or  n  is  often  used  as  the  symbol  for  the  inverse  of 
the  length  scale  factor,  i.e.,  l/Kj. 
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APPENDIX  B:  SCALING  THE  RIGIDICE  MODEL 


The  Rigidice  model  of  frost  heave  is  a  set  of 
differentiEd  equations  expressing  mass,  energy  and 
force  balances  in  the  frozen  fringe  for  air-free,  col¬ 
loid-free  and  solute-free  soil  (Miller  1978).  A  set 
of  reduced  or  scaled  variables  for  this  model  were 
later  presented  by  Miller  (1990).  To  avoid  redun¬ 
dancy,  this  later  work  did  not  contain  derivations 
but  referred  to  the  earlier  classic  works  for  ice- 
free  unsaturated  soil  (Miller  and  Miller  1955,  Miller 
1980a).  For  completeness  of  this  paper,  this  ap¬ 
pendix  presents  a  brief  derivation  of  the  reduced 
variables  presented  by  Miller  for  the  liquid  wa¬ 
ter  and  solid  ice  phases.  It  also  serves  as  an  exam¬ 
ple  of  how  the  seemingly  subjective  process  of 
scaling  equations  is  performed.  All  scaled  vari¬ 
ables  will  be  given  upper-case  letters  to  distin¬ 
guish  them  from  the  nonscaled  dimensional  vari¬ 
ables  that  are  in  lower  case. 


Ice 


Figure  Bl.  Equilibrium 
configuration. 


We  start  with  Young's  equation  relating  the 
contact  angle  a  and  the  surface  tensions  for  all 
three  phases  of  water: 

YiA=Ynv  +  YwACOsa  (Bl) 

in  which  Y  is  the  surface  tension  and  the  subscripts 
I,  W  and  A  are  for  ice,  water  and  air,  respectively. 
Figure  Bl  is  a  schematic  of  the  equilibrium  con¬ 
figuration  and  force  diagram  about  the  point  of 
mutual  contact  for  all  three  interfaces.  Since  cos 
a  is  already  dimensionless,  eq  Bl  can  only  be  made 
dimensionless  by  dividing  through  by  one  of  the 
surface  tensions.  A  reasonable  choice  is  yip^  with 
the  largest  magnitude,  which  results  in  the  re¬ 
duced  Young's  equation 

1  =  Fjyv  +  Fwa  ® 

where  the  reduced  surface  tensions  are 


and 


“  .y  Yaw 
Yia 


(B3) 


To  simplify,  the  remainder  of  this  appendix  will 
address  ice-water  systems.  A  similar  set  of  equa¬ 
tions  for  water-air  and  ice-air  systems  are  ob¬ 
tained  by  direct  analogy  to  the  ice-water  equa¬ 
tion  and  are  listed  by  Miller  (1990). 

The  distinguishing  trait  exhibited  by  water  and 
ice  in  a  porous  medium  that  does  not  exist  in  un¬ 
restricted  bulk  is  the  existence  of  curved  inter¬ 
faces.  The  force  balance  across  a  curved  interface 
separating  two  phases  in  equilibrium  is  given  by 
the  Laplace  surface  tension  equation,  which  for 
an  ice-water  interface  is 


Pi-Pw=2^  (B4) 

where  the  absolute  pressure  is  p  (N/m^)  and  the 
mean  radius  of  curvature  is  r.  In  the  case  of  fro¬ 
zen  soil,  this  equation  expresses  the  behavior  of 
water  and  ice  contained  within  the  pores.  This 
suggests  a  length  scale  on  the  order  of  the  pore 
size,  so  we  introduce  a  microscopic  length  X.  Equa¬ 
tion  B4  contains  three  variables,  pi,  pw  and  r,  and 
one  constant,  Yj^^.  Equation  B3  shows  how  to  re¬ 
duce  the  surface  tension.  At  this  point  we  need 
to  decide  how  to  reduce  the  remaining  variables. 
The  simplest  approach  is  to  acknowledge  that  the 
radius  of  curvature  is  a  pore  length  dimension, 
so  it  is  reduced  by  dividing  by  X.  The  reduced 
curvature  is  therefore 

Substituting  B3  and  B5  into  the  right  hand  side 
of  B4  determines  the  reduced  pressure  by  dimen¬ 
sional  equivalence: 

P^—p  (B6) 

Yia 

giving  the  reduced  form  of  the  Laplace  surface 
tension  equation: 

Kiw 

At  thermodynamic  equilibrium  the  chemical 
potentials  of  ice  and  water  are  equal  and  are  ex- 
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pressed  by  the  differential  form  of  the  Gibb's  re¬ 
lationship.  It  has  been  shown  (Loch  1978,  Black 
1992)  that  the  Gibb's  relationship  for  monocom¬ 
ponent  and  multicomponent  systems  contained 
within  a  porous  medium  that  subjects  the  water 
to  a  general  body  force  can  be  expressed  in  a  sim¬ 
plified  Clapeyron  equation.  When  pressures  are 
given  in  standard  gauge  pressures  w,  the  Clapey¬ 
ron  equation  is  written  as 

=  (B8) 

IlW 

in  which  hiyj  is  the  volumetric  latent  heat  of  fu¬ 
sion  (J/m^),yjy^  is  the  specific  gravity  of  ice  and 
6  (°C)  is  the  elevation  of  the  melting  point  of  ice 
from  the  standard  melting  point  of  ice  Bq.  For 
normal  water,  0  will  be  negative. 

The  gauge  pressures  will  scale  the  same  as  the 
absolute  pressures: 

Us—u.  (B9) 

Yia 

The  specific  gravity  of  ice  is  already  dimension¬ 
less,  so  the  remaining  constants  Gq  and  hiyj  must 
be  used  to  reduced  the  only  variable:  tempera¬ 
ture  0.  Earlier  scaling  attempts  by  Miller  (Black 
1985)  developed  a  reduced  temperature  that  in¬ 
cluded  the  latent  heat  term,  but  for  the  benefit  of 
clarity.  Miller  (1990)  suggested  a  different  reduced 
temperature: 

©a±0.  (BIO) 

Substituting  these  scaling  parameters  for  pres¬ 
sure  and  temperature  into  eq  B8  results  in  an  ad¬ 
ditional  reduced  term  called  the  reduced  latent 
heat: 

(Bll) 

Yia 

These  scaling  relationships  give  a  reduced 
Clapeyron  equation: 


The  Laplace  and  Clapeyron  equations  show 
that  the  pressures,  temperature  and  curvature  are 
interrelated.  When  the  porous  medium  is  geo¬ 
metrically  complicated,  as  in  the  case  of  soil,  it  is 
helpful  to  simplify  the  pressures-curvature  be¬ 
havior  by  expressing  the  Laplace  equation  with 
empirically  determined  material  properties.  In  our 
case  the  necessary  material  properties  are  the  vol¬ 
ume  fractions  for  water  W{uy/,  Wj),  ice  Uj) 


and  grains  G,  These  volume  fractions  are  related 
by 

G4*  W(w^,Ui)  + =  1.  (B13) 

The  complete  scaled  form  is  simply 

G  +  W(L[w ,  )  +  l(Ziw .  )  =  1  •  (B14) 

The  criterion  for  lens  initiation  is  obtained  from 
the  geotechnical  concept  of  effective  stress,  or  the 
particle-to-particle  contact  stress.  When  this  stress 
is  zero,  the  particles  no  longer  remain  in  contact 
and  an  ice  lens  is  able  to  form.  The  Snyder-Terzaghi 
equation  partitions  the  total  confining  stress  Oj 
(N/m^)  between  the  effective  (N/m^)  and  neu¬ 
tral  or  pore  (N/m^)  stresses  at  locations  near  the 
initiation  of  an  ice  lens: 

aT=Oe+an.  (B15) 

Figure  B2  shows  the  partitioning  of  these  stresses. 
The  effective  and  neutral  stresses  in  this  relation¬ 
ship  are  different  from  the  traditional  stresses  in 
the  Terzaghi  equation  because  these  include 
Snyder's  flawed  solid  theory  (Snyder  and  Miller 
1985).  These  stresses,  like  the  pressure,  take  place 
at  the  pore  scale,  so  they  will  scale  just  as  pres¬ 
sures  do: 

^  X 

1  = - a.  (B16) 

Yia 

The  flow  of  water  through  the  frozen  fringe  is 
assumed  to  obey  Darcy's  law: 

i^vv  =  (^w/^i)(/w  (B17) 

in  which  Vy^  is  the  volumetric  flux  of  water 
(m^/ m^s),  Uj)  is  the  capillary  conductivity 

(m^/Ns),  fyj  is  the  body  force  per  unit  volume 
(N/m^)  and  V  is  the  gradient.  To  reduce  Darcy's 
law  we  begin  by  examining  the  last  term.  This  is 
the  gradient  of  water  pressure  across  the  region 
of  interest.  If  we  are  looking  at  flow  through  an 
individual  pore,  then  this  length  would  appear 
to  be  X.  This  would  be  correct  for  this  special  and 
very  limited  case.  In  general,  far  more  complex 
flow  patterns  through  a  myriad  of  pore  geom¬ 
etries  containing  different  amounts  of  ice  and 


Figure  B2.  Stress  partition¬ 
ing  of  ice,  water  and  air. 
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water  are  of  concern.  At  this  scale  we  refer  to  the 
porous  medium  behaving  as  a  continuum  display¬ 
ing  rmiform  material  properties,  and  we  are  inter¬ 
ested  in  the  flow  into  and  then  out  of  a  volume 
of  material.  This  suggests  that  a  different  length 
scale  is  required  to  describe  the  gradient,  since  X 
is  applicable  only  at  the  pore  scale.  The  additional 
macroscopic  length  scale  is  Using  this  new 
length  scale  to  reduce  the  gradient  operator  re¬ 
sults  in  the  reduced  body  force: 


P  -  f 

%  — /w 

'AW 


(B18) 


The  proportionality  factor  in  Darcy's  equation 
is  the  source  of  much  confusion.  In  eq  B17  it  is 
expressed  as  capillary  conductivity  Soil  sci¬ 
entists  tend  to  use  hydraulic  conductivity  K  (m/ 
s),  and  geologists  and  petroleum  engineers  use 
intrinsic  permeability  k  (m^).  These  three  param¬ 
eters  are  related  by 


k 

h 


(B19) 


in  which  p  (kg/m^)  is  the  density  of  water,  g  is 
the  acceleration  due  to  gravity  (or  other  body 
force)  and  Ti  (N/m^s)  is  the  viscosity  of  water. 
They  are  all  functions  of  water  and  ice  content, 
which  are  functions  of  the  water  and  ice  pres¬ 
sures. 

The  process  of  scaling  the  capillary  conduc¬ 
tivity  consists  of  first  realizing  that  it  is  already 
inversely  proportional  to  viscosity  (i.e.,  eq  B19), 
so  it  must  therefore  be  scaled  by  multiplying  by 
the  viscosity.  It  then  becomes  the  intrinsic  perme¬ 
ability,  which  is  in  dimensions  of  length  squared. 
To  complete  the  scaling  we  must  decide  on  the 
correct  length  scale  to  use.  Again  relying  on  the 
continuum  approximation,  the  conductivity 
should  not  depend  on  the  size  of  the  sample  once 
it  attains  the  minimum  size  required  for  the  con¬ 
tinuum  approximation.  This  requires  the  use  of 
the  microscopic  length  scale  X  to  scale  the  con¬ 
ductivity.  Miller  and  Miller  (1955)  deduced  the 
same  length  scale  choice  by  considering  the 
Navier— Stokes  equation  applied  to  flow  through 
the  pore.  The  reduced  capillary  conductivity  is 
then 


(B20) 

Substituting  eq  B20,  B19,  B6  and  the  scaling  fac¬ 
tor  ^  to  reduce  the  gradient  operator  into  eq  B17 
gives  the  scaled  form  of  Darcy's  equation: 


(B21) 

in  which  the  reduced  volumetric  flux  of  water 
Vyf  is,  by  balancing  constants, 

V'w-A-w  (B22) 

Mia 

The  continuity  equation  for  pore  constituents 
is 


V  •  (i7^  + 

=  -  •”  [  (ww  /  )  +  Yiw  f  (^w  /  )  ]  (B23) 

in  which  Vi  (m^/m^s)  is  the  volumetric  flux  of 
ice  and  t  is  time.  The  operator  V«  is  the  diver¬ 
gence,  which  simply  scales  by  multiplying  by 
The  volumetric  flux  of  ice  must  scale  like  the 
volumetric  flux  of  water  to  maintain  dimensional 
equality  on  the  left  side  of  eq  B23.  Since  the  volu¬ 
metric  water  and  ice  contents  are  already  dimen¬ 
sionless,  the  reduced  time  T  is 

(B24) 

The  conservation  of  thermal  energy  is 

+  liiwW(uw/Mi)]  (^25) 

in  which  uh  (W/m^)  is  the  macroscopic  flux  of  sen¬ 
sible  heat  and  Ch(mw/  “r)  K)  is  the  volumet¬ 
ric  heat  capacity  excluding  heat  resulting  from 
phase  change.  Tlie  volumetric  heat  capacity  is  as¬ 
sumed  to  be  a  material  property  that  is  constant 
in  the  continuum  assumption  for  a  given  state  of 
ice  and  water  contents.  Starting  with  the  left  side 
of  eq  B25  the  reduced  macroscopic  flux  of  sen¬ 
sible  heat  l/jj  is  obtained 

yjA 

The  reduced  volumetric  heat  capacity  must 
therefore  be 


Cr  (^w  /  ^i)  =  (^w  /  ^i)  •  (B27) 

flA 

The  macroscopic  flux  of  sensible  heat  is  given 
by  Fourier's  law: 
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(^W  /  )  V0  (B28) 

in  which  kn  (W/mK)  is  the  thermal  conductiv¬ 
ity.  This  conductivity  is  also  a  material  property 
that  is  assumed  to  be  constant  in  the  continuum 
sense  for  a  given  state  of  water  and  ice  contents. 
The  reduced  thermal  conductivity  is  then 

Kh  {Uyj,Ui)  =  [uy^,Ui).  (B29) 

yiA 

The  most  important  characteristic  of  the  Rigidice 
model  that  sets  it  apart  from  others  is  the  concept 
of  thermally  induced  regelation.  This  results  in  the 
inclusion  of  an  extra  mass  flux  resulting  from  the 
bulk  movement  of  ice  in  the  frozen  fringe.  The  ve¬ 
locity  of  ice  movement  in  the  fringe  is  the  heave 
velocity  (m/s)  relative  to  the  stationary  grains 
of  the  porous  media  below  the  fringe.  The  volu¬ 
metric  flux  of  ice  Vi  (m^/m^s)  resulting  from  this 
regelation  is 

Vi  =  (^30) 

From  eq  B22  we  know  how  volumetric  fluxes  scale, 
so  the  reduced  regelation  velocity  is  therefore 

(B31) 

MIA 

This  completes  the  scaling  of  the  more  impor¬ 
tant  variables  and  parameters  involved  in  ground 
freezing.  A  list  of  all  scaled  variables  applicable 
to  ground  freezing,  as  well  as  all  water  phases, 
was  presented  by  Miller  (1990). 

An  example  of  how  heaving  experiments  should 
be  interpreted  by  reduced  variables  seems  in  or¬ 
der.  Given  two  experiments  using  the  same  soil 


but  subjected  to  different  body  forces,  eq  B18,  BIO 
and  B24  show  how  they  must  relate.  Using  the 
convention  of  calling  the  laboratory  reference 
experiment  the  prototype  p  and  the  experiment 
subjected  to  the  different  body  force  m  for  model, 
the  reduced  body  force  f ^  must  be  equal  for  both: 

(B32) 

"lA  ^lA 

If  we  say  that/j^  =  100 /p,  as  might  be  expected  in 
a  centrifuge,  and  we  use  the  same  soil  in  both 
experiments,  then  the  total  lengths  of  the  two  ex¬ 
periments  are  related  by  =  0-01  Cp-  Since  the 
reduced  temperatures  0  must  be  equal  at  corre¬ 
sponding  locations,  the  actual  temperatures  0  will 
therefore  be  the  same  in  both  the  prototype  and 
the  model.  This  leads  to  the  most  important  ex¬ 
perimental  consequence  of  centrifuge  modeling. 
Since  the  reduced  times  T  must  be  equal,  the  ac¬ 
tual  times  are  related  by  =  t^/100^.  Scaled  ex¬ 
periments  are  finished  in  a  very  timely  manner 
compared  to  similar  laboratory  experiments. 

A  final  note  is  on  the  meaning  of  the  micro¬ 
scopic  length  X.  If  the  porous  medium  consisted 
of  identical  pores,  then  this  length  would  just  be 
the  pore  size.  But  in  real  porous  media  there  is  a 
distribution  of  sizes  and  geometries  of  pores.  The 
microscopic  length  in  these  cases  represents  some 
averaged  value  of  the  pores.  It  is  more  like  the 
size  of  the  smallest  representative  elementary 
volume  (REV)  that  begins  to  describe  the  porous 
media  as  a  continuum  of  uniform  material  prop¬ 
erties.  The  difficulty  is  in  actually  measuring  this 
parameter  without  completely  destroying  the 
original  material. 
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